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Abstract Engineered wood flooring (EWF) is a multi¬ 
layer composite flooring product. The cross layered struc¬ 
ture is designed to give good dimensional stability to the 
EWF under changing environmental conditions. However, 
during winter season in North America, the indoor relative 
humidity could decrease dramatically and generate an 
important cupping deformation. The main objective of this 
study was to characterize the interlaminar stresses (cr 33 , cr 13 
and ( 723 ) distribution at free-edges in EWF made with an 
OSB substrate. A three-dimensional (3D) finite element 
model was used to predict the cupping deformation and to 
characterize stresses developed in the EWF. The finite 
element model is based on an unsteady-state moisture 
transfer equation, a mechanical equilibrium equation and 
an elastic constitutive law. The physical and mechanical 
properties of OSB substrate were experimentally deter¬ 
mined as a function of the density and moisture content. 
The simulated EWF deformations were compared against 
the laboratory observations. For both simulation and 
experimental results, the cupping deformation of EWF was 
induced by varying the ambient relative humidity from 50 
to 20% at 20°C. A good agreement has been found between 
the numerical and experimental EWF cupping deformation. 
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The stress distribution fields generated by the model cor¬ 
respond to the delaminations observed on the OSB sub¬ 
strate in the climate room. Delamination in EWF can occur 
principally under the action of the tension stress or a 
combination of tension and shear stresses. Finally, simu¬ 
lated results show that the levels of interlaminar stresses 
are maximal near the free-edges of EWF strips. 

Keywords Layered structure • Delamination • 

Finite element analysis (FEA) • Wood testing 


Introduction 

Engineered wood flooring (EWF) is a cross layered structure 
designed to give good dimensional stability under changing 
environmental conditions. However, during winter season in 
North America, the indoor relative humidity could decrease 
dramatically and generate an important cupping deforma¬ 
tion. The magnitude of EWF deformation depends on the 
physical and mechanical properties of each layer and its 
behavior following changes in moisture content (MC) and 
temperature near the top surface. 

Barbuta et al. [1] developed two types of OSB (Oriented 
Strand Board) as a substitute for Baltic Birch (Betula 
pendula Roth) plywood (BBP) used as a substrate for 
EWF. The specialty OSB panels were manufactured from 
two types of wood strands: a mixture of 90% aspen 
(Populus tremuloides Michx.) and 10% paper birch (Betula 
papyrifera Marsh.) and 100% ponderosa pine (Pinus pon- 
derosa Dougl. ex Laws.). In a second study, Barbuta et al. 
[2] evaluated the performance of EWF prototypes manu¬ 
factured with OSB substrate in conditioning room. They 
found that there was no significant difference, in terms of 
cupping deformation, between the EWF manufactured with 
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aspen/birch OSB and BBP substrate. However, after the 
evaluation test in conditioning room, delaminations were 
observed in the specialty OSB substrate. These delamina¬ 
tions are the result of the high level of interlaminar stresses 
distribution at free-edges in EWF. 

The finite element method is a useful tool to predict the 
behavior of EWF [3, 4]. Blanchet et al. [5] developed a 
three-dimensional (3D) finite element model of the 
hygromechanical cupping in layered wood composite floor¬ 
ing. The model was based on two sets of equations: (1) the 
three-dimensional equations of unsteady-state moisture 
diffusion and (2) the three-dimensional equations of elas¬ 
ticity including the orthotropic Hooke’s law which takes 
into account the shrinkage and swelling of each layer. 
Blanchet et al. [6] used this finite element model to support 
the design of EWF. They first demonstrated that at the 
macroscopic level, the hypothesis of elastic linearity is 
acceptable. Geometrical parameters and wood species 
properties were then used to assess different EWF con¬ 
struction. They found that the thickness and mechanical 
properties of the core layer are the most important design 
parameters controlling cupping deformation of EWF in 
service. They also indicated that the backing layer has a 
smaller, although still significant, impact on cupping. 

The wood-adhesive interface was characterized by 
Belleville et al. [7] to determine its impact on EWF hy¬ 
gromechanical behavior. The finite element analysis has 
shown that better results are obtained when the wood- 
adhesive interface is considered in the model. 

Deteix et al. [8] have investigated the impact of the 
adhesive layer on the dimensional stability of layered wood 
composites and how it should be modeled by the finite 
element method. They showed that the use of a finite ele¬ 
ment mesh with a single layer of elements in the adhesive 
layer can lead to significant errors. This is particularly true 
for low values of the adhesive effective diffusion 
coefficient. 

The objectives of the current study are to determine the 
physical and mechanical properties of the OSB substrate 
and to characterize stresses leading to delamination at the 
free-edge of EWF with an OSB substrate by the finite 
element method. 


Materials and methods 

Physical and mechanical properties of EWF 
components 

The mechanical and physical properties of OSB panels 
were experimentally determined in laboratory. Unidirec¬ 
tional OSB panels were manufactured with three target 

Q 

densities: 500, 650 and 800 kg/m . Flat density profiles 


through panel thickness were achieved using a temperature 
of 100°C during press closing. When the nominal thickness 
was reached, the platens temperature was raised until the 
core temperature achieved 140°C. The platen temperature 
corresponding to a core temperature of 140°C was 165°C. 
The other OSB manufacturing parameters were the same as 
those used by Barbuta et al. [1, 2] in the production of 
aspen/birch OSB panels. A total of 30 panels were manu¬ 
factured. The physical and mechanical properties were 
determined for each density. The mechanical properties 
were also determined for two RH levels: 50 and 20% at 
20°C. These two environmental conditions represent the 
initial and final stages of the numerical simulation. The 
properties determined were modulus of elasticity (F L , E T 
and F r ), shear modulus (G LR and G TR ), interlaminar shear 
strength (S RR and S TR ), internal bond strength (IB), 
shrinkage coefficients (/? L , /? T and /? R ) and the diffusion 
coefficients (D L , D T and D R ). Tension and shear properties 
were determined using the ASTM D 1037 9 standard. The 
F, T and R directions correspond to the longitudinal, 
transverse and through-thickness directions of the unidi¬ 
rectional OSB panel. The modulus of elasticity perpen¬ 
dicular to the panel surface (F R ) was determined using a 
5 mm MTS extensometer (Model 632.29F-30) attached to 
the test specimen. The shrinkage coefficients for each 
density between 50% RH and 20% RH were determined 
using Eq. 1: 

_ d 5 Q — d 2 o 
dsoAM 

where d 2 o is the dimension at 20% RH (mm), d 50 is the 
dimension at 50% RH (mm), and AM is the moisture 
content difference (%) 

The steady-state method used to determine the diffusion 
coefficient was inspired from Siau [11]. Twenty-five dif¬ 
fusion cells were prepared for each density. The water 
vapor diffusion coefficients were determined by the fol¬ 
lowing equation: 

lOOwL N 

D =-— (2) 

where w/t is the water vapor flow (kg s -1 ), L is the length 
in the flow direction (m), A is the cross-section area of 
specimen normal to the flow direction (nT), p w is the water 

_o 

density (1000 kg m ), G is the specific gravity of wood at 
the average moisture content, and AM is the moisture 
content difference across specimen thickness (%) 

To model the behavior of EWF, the physical and 
mechanical properties of all materials were needed. The 
properties of maple sugar layer, PVA adhesive type I and 
the six Poisson’s ratios (v;y) for OSB were also taken from 
the literature. The parameters used in the numerical model 
are presented in Table 1. 
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Finite element modeling 


The EWF strip considered in this study is made of a 3-mm 
thick sugar maple surface layer and a 9-mm thick OSB 
substrate layer (Fig. 1). The OSB substrate is composed of 
three cross-oriented layers. The ratio between the two 
surface layers and the core layer is 0.45/0.1/0.45. The 
adhesive used to bond the surface layer and the substrate is 
a type I PVA (polyvinyl acetate). The thickness of the 
adhesive layer was estimated to be 0.1 mm [8]. The 
construction of EWF strips was designed such that the longi¬ 
tudinal axis of the substrate surface layers would be 
perpendicular to that of the sugar maple wood layer. 

The cupping deformation was generated by varying 
relative humidity (RH) from 50 to 20% at 20°C. This 
corresponds to a decrease in moisture content from 8.6 to 
5.8% for the sugar maple surface layer and PVA adhesive 
and from 6.2 to 3.4% for OSB substrate. The water vapor 
desorption occurred by free convection at the top of the 
surface layer. All the other surfaces were assumed 
impervious. 

The theoretical model of hygromechanical EWF defor¬ 
mation used in this study is based on the model developed 
by Blanchet et al. [5] and improved by Deteix et al. [8] and 
Belleville et al. [7]. The governing equations used in the 
model are three-dimensional moisture conservation (3 and 
4), three-dimensional equation of mechanical equilibrium 
(5) and Hooke’s law for orthotropic materials which takes 
into account the shrinkage of each layer (6). 

The moisture transfer is described by the equation: 


db 0M 

Ioo"a7 


+ v- 


(-4vm) 



with K m = 


Dd \, 

Too 



where K M is the tensor of effective water conductivity 
(kg m -1 s' 1 % _1 ), D : tensor of effective moisture diffu- 

o _i _o 

sion (nT s ), d h is the basic density (kg m ), t is the time 
(s) and M is the moisture content (%). 

The moisture transfer in EWF can be written in the 
following manner: 


0M J 0 0 0 
dt \ dx\ 0X2 0X3 

( 4 ) 




= 0 


The mechanical model is based on the three-dimensional 
equation of mechanical equilibrium: 

g = 0 i,7 =1,2,3 (5) 

where is the normal and shear stress components (Pa). 

It is assumed that the wood surface layer, the OSB 
substrate and the adhesive are elastic. The stress tensor is 


calculated using Hooke’s law for the stress-strain 
relationship. 

®ij Cijkl(£kl ( 6 ) 

where C ijk i is the stiffness tensor, 8 ki is the normal and shear 
strain components, p k i is the moisture shrinkage coeffi¬ 
cients (% _1 ), and AM is the moisture content change 
between two time steps (%). 

The normal and shear strains and displacements u t are 
related by the equation: 

1 (dui duj\ 

£ij = 2 ( v 0^ + 0^j 

The finite element (FE) modeling of hygromechanical 
deformation of EWF was performed using the FE code 
MEF++ developed by the GIREF (Interdisciplinary 
Research Group on Finite Elements) at Universite Laval, 
Quebec, Canada. The computational domain is a 50-mm 
length EWF strip (Fig. 2). This dimension corresponds to 
EWF sample length generally used in the delamination test. 
Two FE meshes corresponding to half of this domain were 
generated with increasing element densities. The relatively 
coarse and fine meshes contained 117900 and 495450 
ten-node tetrahedral elements. The 0.1-mm thick adhesive 
layer was meshed with three layers of elements. Figures 3 
and 4 illustrate the coarse and the fine meshes used in this 
numerical experiment. 

Each layer of EWF is assumed to be elastic. It is also 
assumed that the principal material directions of the sugar 
maple surface layer and the three sublayers of the OSB 
substrate are perfectly oriented with the strip. The experi¬ 
mental data obtained for unidirectional OSB panels were 
linearly extrapolated to calculate the physical and 
mechanical properties of OSB substrate function of the 
vertical density profile (VDP) and the AM. For FE simu¬ 
lations, the initial moisture content M 0 was set at 8.6% for 
the sugar maple surface layer and PVA adhesive and 6.2% 
for OSB substrate. This was due to the different sorption 
isotherm of these two materials. The equilibrium moisture 
content was as follows: 5.8% for surface layer and 
PVA adhesive and 3.4% for OSB substrate. It was assumed 
that moisture transfer occurs only through the top of the 
surface layer. 

Boundary condition: 

q = h(M — Mqo) on exposed face (8) 

q = 0 on the other surfaces (9) 

_o _ 1 

where, q is the moisture flux (kg water m wet wood s ), h is the 

_o _i _i 

convective mass transfer coefficient (kg m s % ), 

M is the moisture content (%), and is the moisture 
content at equilibrium (%) 

The initial condition for moisture content was: 



4T Springer 



















330 


J Wood Sci (2012) 58:327-335 


Table 1 Finite element model parameters 

Parameters 

Material 




















Surface 

Substrate 
















Adhesive 


Sugar maple 

OSB 

















PYA I 

d h (kg/m 3 ) 

597 a 


















1178 1 

d„ (kg/m 3 ) 

D (m 2 s -1 ) 


500 f 





650 f 






800* 






4.18 x 10“ 12i 

D l (m 2 s- 1 ) 

2.20 x 10 _9b 

— 





— 






— 







D t (m 2 s _1 ) 

1.80 x 10~ llc 

— 





— 






— 







D r (m 2 s- 1 ) 

1.80 x 10“ llc 

2.63 x 

10* 11 




1.33 

X 

10-“ f 




6.68 

X 

10 _12f 





r(°c) 


20 f 





20 






20 







RH (%) 


50 


20 



50 



20 



50 



20 




D r (m 2 s' 1 ) 

1.80 x 10 _llc 

2.63 x 

10 _llf 




1.33 

X 

10-“ f 




6.68 

X 

10 _12f 





M 0 (%) 

8.6 

6.2 


— 



6.2 



— 



6.2 



— 



8.6 

M x (%) 

5.8 

— 


3.4 



— 



3.4 



— 



3.4 



5.8 

h (kg m -2 s - 1 
% _1 ) 

P (m m _1 % _1 ) 

3.20 x 10 _4b 






3.20 

X 

10- 4b 










3.20 x 10 _4b 
3 x 10“ 3i 

ft. (m m _1 
% _1 ) 

1.50 x 10 _4d 

2.36 x 

10 _4f 




2.07 

X 

10~ 4f 




1.75 

X 

10 _4f 





ft- (m m“ 

3.30 x 10 _3d 

9.18 x 

l 0 - 4f 




9.14 

X 

10~ 4f 




9.79 

X 

i0 - 4f 





Pr (m m -1 
% _1 ) 

2.10 x 10 _3d 

2.86 x 

10- 3f 




3.57 

X 

10- 3f 




5.00 

X 

io- 3f 




1.27 x 10 +loj 

E (Pa) 

E R (Pa) 

1.38 x 10 +1Oe 

6.72 x 

10 +9f 

7.00 

X 

10 +9f 

9.50 

X 

10 +9f 

9.83 

X 

10 +9f 

1.23 

X 

10 +lof 

1.30 

X 

10+i° f 


E t (Pa) 

6.78 x 10 +8e 

9.27 x 

10 +9f 

9.84 

X 

10+8f 

1.54 

X 

10 +9f 

1.61 

X 

10 +9f 

2.15 

X 

10 +9f 

2.21 

X 

10 +9f 


E r (Pa) G (Pa) 

1.31 x 10 +9e 

7.50 x 

10 +7f 

8.50 

X 

10 +7f 

1.26 

X 

10+8f 

1.51 

X 

l 0 + 8f 

1.47 

X 

10+8 f 

1.70 

X 

l 0 + 8f 

4.72 x 10 +9k 

G lr (Pa) 

1.01 x 10 +9e 

1.54 x 

10+8 f 

1.62 

X 

10 +8f 

2.24 

X 

10+8 f 

2.29 

X 

10+ 8f 

2.45 

X 

10+8 f 

2.44 

X 

10+ 8f 


G T r (Pa) 

2.55 x 10 +8e 

7.70 x 

10 +7f 

8.20 

X 

10 +7f 

1.06 

X 

10+8f 

1.36 

X 

io+s f 

1.24 

X 

10+Sf 

1.47 

X 

10+8f 


Glt (Pa) v 

7.53 x 10 +8e 

3.37 x 

10 +9g 

1.12 

X 

10 +9g 

1.63 

X 

10+ 9g 

1.99 

X 

l 0 + 9g 

2.57 

X 

10 +9g 

2.68 

X 

l 0 + 9g 

0.35 1 

V LT 

0.500 e 

0.183 h 


















Vrt 

0.820 e 

0.019 h 


















V TL 

0.025 e 

0.161 h 


















Vrl 

0.044 e 

0.013 h 


















V T r 

0.420 e 

0.312 h 


















V LR 

0.460 e 

0.364 h 



















d h basic density, d n nominal density, D z diffusion coefficient in the i direction, T temperature, RH relative humidity, M 0 initial moisture content, M\ 
final moisture content, h convective mass transfer coefficient, ft shrinkage coefficient in the i direction; ft elastic modulus in the i direction; Gy shear 
modulus in the ij plane; Vy Poisson’s ratio in the ij plane 

a Jessome [10] 

b Siau [11] 

c Blanchet et al. [5] 


Goulet and Fortin [12] 


e Bodig and Jayne [13] 

-C 

Experimental value 
g Estimated value 
h Zhu et al. [14] 

1 Belleville et al. [7] 

J Bandrup et al. [15] 
k Ballon and Drolet [16] 
1 Urayama et al. [17] 
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Fig. 1 Construction and 
geometry of EWF prototype 


Fig. 2 Domain geometry and 
boundary conditions 





Fig. 3 Coarse finite element 
mesh used in the numerical 
model 



Fig. 4 Fine finite element mesh 
used in the numerical model 
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• for the surface and the adhesive layers 

M(xi,* 2 ,* 3 , to) = M 0 = 8 . 6 % V(xi,x 2 ,x 3 ) 

( 10 ) 

• for the OSB substrate 


M(xi,X 2 ,x 3 ,to) = M 0 = 6.2% V(xi,x 2 ,x 3 ) 

( 11 ) 

The moisture content at equilibrium was: 

• for the surface and the adhesive layers 

M(x 1 ,x 2 ,x 3 ,t 00 ) =M ao =5.8% \/{x\,x 2 ,x 3 ) 

( 12 ) 

• for the OSB substrate 

M(x u x 2 ,x 3 ,t 00 ) = Moo = 3.4% V(xi,x 2 ,x 3 ) 

(13) 


For the mechanical part of the problem, the conditions 


were: 

Initial condition: 

o(xi,X2,X'i,to) = (To = 0 V(xi,x 2 ,x 3 ) (14) 

u(xi,X 2 ,X 3 ,t 0 ) = Mo = o V(x l,x 2 ,x 3 ) (15) 

Boundary conditions (Fig. 2): 
u\ = 0 at ( 0 ,^ 2 , x 3 ) (16) 

u 2 = 0 at (0,43,0) and (0,43,0.675) (17) 

u 3 = 0 at (0,43,0) (18) 


The numerical simulations were performed using a 
constant time step of 0.1 s. The total time of simulation 
was 28 days. At each time step, the moisture content 
variation (AM) was used to calculate the corresponding 
displacements ( u ). The AM was also used to update the 
mechanical and physical properties of OSB substrate. To 
reduce the computation time and the data generated by the 
model, the mechanical calculations were realized at only 
20 time steps. 

To compare the simulated EWF cupping deformations 
with the experimental deformations, EWF prototypes were 
produced. The manufacturing parameters are the same as 
those used by Barbuta et al. [2] for EWF prototypes with an 
OSB substrate. After a conditioning period at 50% RH and 
20°C, twelve EWF specimens of 86 mm wide and 50 mm 
long were placed in a conditioning room at 20% RH and 
20°C for 4 weeks. Before the test in climate chamber, the 
EWF specimens were sealed with aluminum foil on the 
bottom face and the edges to impose moisture transfer by 
the top surface only. The cupping deformations were 
measured over the width of the strip with a dial gauge as 
described by Blanchet et al. [18]. The coordinates of 
measurement points (A, B , C) were introduced in the finite 
element model to calculate the simulated cupping defor¬ 
mations (Fig. 2). The cupping deformation was defined as 
the difference in vertical displacement between the points 
A and B and the central point C. The minus sign was used 
to indicate that the cupping is a convex deformation of the 
EWF strip. 


Results and discussion 

Engineered wood flooring cupping deformation 

Figure 5 presents experimental and numerical cupping 
deformation as a function of time. The simulated curves 
from the coarse and fine meshes show a good agreement 
with experimental data. We note that the coarse mesh 
provides results closer to the experimental measurements. 
However, the difference between the results generated by 
the two meshes is small. Compared to the experimental 
data, the fine mesh solution overestimated the maximum 
deformation by 7%. 

Even though the cupping simulated curve obtained from 
a coarse mesh is closer to the experimental one, a finer 
mesh is necessary for an accurate description of the inter¬ 
laminar normal and shear stresses. For both experimental 
and numerical results (Fig. 5), approximately 91% of the 
maximum cupping deformation occurred in the first 
10 days of the experiment. The slopes of the simulated 
curves in the first days are influenced by the properties of 
the surface layer. Blanchet et al. [5] showed that the vari¬ 
ation of the transverse water diffusion coefficient has an 
important effect on the curve slope. After that period, the 
cupping deformation curve turns into a relatively flat pla¬ 
teau. It is important to indicate that the numerical simula¬ 
tions were performed using the mean values of the material 
properties without considering their variability. The aver¬ 
age coefficients of variation (CV) for physical and 
mechanical properties of clear wood are between 15 and 
25% [19], whereas the average measured CV values for 
OSB are approximately 30%. Parametric studies indicate 
that the maximum cupping deformation is influenced by 
the modulus of elasticity in parallel and perpendicular to 
the surface grain of OSB and by sugar maple shrinkage 



Fig. 5 Cupping deformation as function of time. Experimental and 
numerical results 
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Fig. 6 Comparison between the calculated normal stress cr 33 at the 
tongue and groove edges after 28 days at 20°C and 20% RH and the 
internal bond strength IB of the OSB substrate at 20% RH as function 
of 

coefficient in tangential direction. The variation of these 
parameters with ±15% leads to a modification of the 
simulated maximum cupping deformation of approxi¬ 
mately ±19%. 

Stress characterization 

Near the free-edges of the EWF strip, interlaminar normal 
and shear stresses rise sharply. These interlaminar stresses 
are mainly responsible for delamination fracture. In our 
case, the interlaminar stresses were generated by the mis¬ 
match in the shrinkage coefficients and elastic properties 
between the surface layer and the substrate during moisture 
movement. Because the shrinkage coefficient of the sugar 
maple surface layer is very small in the longitudinal 
direction (Table 1), a low level of shear stress cr 13 is 
expected. The maximum of the simulated shear stress was 
0.01 MPa. The distributions of interlaminar normal er 33 and 
shear er 23 stresses at the tongue and groove edges at 
Xi = 12 mm, function of x 3 , are presented in Figs. 6 and 7. 

The two curves in each figure represent the calculated 
values of the stresses at the end of the 28 days experiment 
at 20°C and 20% RH. As we can see, the stresses cr 33 and 
g 2 3 rise significantly in the OSB substrate at the vicinity of 
the adhesive line (v 3 = 9 mm). The peak normal stress er 33 
on the tongue side is higher than that on the groove side 
[2.88 MPa compared to 2.11 MPa (Figs. 6, 8)]. This dif¬ 
ference is attributed to the presence of the groove in the 
OSB substrate. On the groove side, another peak stress is 
generated at v 3 = 7 mm which decreases slowly to the 
bottom of the EWF strip. On the other side, the level of 
maximum stress remains relatively constant until the 



Fig. 7 Comparison between the calculated shear stress cr 23 at the 
tongue and groove edges after 28 days at 20°C and 20% RH and shear 
strength Srr of the OSB substrate at 20% RH as function of x 3 

tongue and decreases to practically zero at v 3 = 7 mm. The 
stress peaks cr 23 on the tongue and groove sides are 
opposite and relatively symmetrical with a maximum of 
2.88 MPa (Figs. 7 and 8). The shear stresses decreased 
rapidly after only 0.3 mm. As in the case of the stress <t 33 , 
the presence of the groove generates another stress peak at 
x 3 = 1 mm. A rough comparison between the simulated 
maximum stresses and the ultimate strength determined 
according to ASTM D 1037 [9] were made. The mean 
values of internal bond strength (IB) and shear strength 
(Slr) as a function of density for unidirectional OSB panels 
with flat density profile (Table 2) were linearly extrapo¬ 
lated to determine the shear and internal bond strengths as a 
function of vertical density profile of the OSB substrate 
(Figs. 6, 7). 

The comparison shows that only the simulated tensile 
stress cr 33 is higher than the internal bond strength 
(Fig. 6). It seems that the model overestimates the stress 
field in the vicinity of the EWF strip free-edges. The use 
of an elastic model may explain this overestimation of 
the stress. Moses and Prion [20] also used an elastic 
finite element model to predict the strength of laminated 
veneer lumber (FVF). They found that shear stresses 
obtained from the model were approximately 1.7 times 
higher than nominal values. 

The stress fields (Fig. 8) obtained from the model are 
in accordance with the observation made on the OSB 
samples after conditioning at 20°C and 20% RH. Seven 
out of twelve samples presented micro-delamination: five 
between the glue line and the tongue and six in the groove 
zone. 

Figure 9 presents the distribution of the calculated normal 
and shear stresses at free-edges function of v 2 . The maximum 
values at free-edges were found at v 3 = 8.34 mm on the 
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Fig. 8 Stress fields calculated by the model; a groove side normal stress cr 33 ; b tongue side normal stress cr 33 ; c groove side shear stress cr 23 ; 
d tongue side shear stress cr 23 


Table 2 Internal bond strength and shear strength of unidirectional 
OSB with flat density profile as a function of nominal density at 20°C 
and 20% RH 


Parameters 

Unidirectional OSB nominal density 



500 kg/m 3 

650 kg/m 3 

800 kg/m 3 

IB (MPa) 

0.583 

0.829 

1.054 

S LR (MPa) 

2.714 

3.898 

4.708 


IB internal bond strength, y LR shear strength in the LR plane 


tongue side and at v 3 = 8.5 mm on the groove side for the 
normal stress 033 and at v 3 = 8.8 mm for shear stress <723 on 
both sides (Fig. 8 ). As expected, the magnitude of the 
interlaminar normal and shear stresses increases rapidly as 
the free-edges are approached. This phenomenon is more 
evident in the case of the normal stress which reaches the 
maximum in approximately 0.5 mm from the free-edge. The 
shear stress begins to increase/decrease at about 30 mm from 
the free-edges. The positive or negative signs of the shear 



Fig. 9 Distribution of the maximum calculated normal cr 33 and shear 
cr 23 stresses at free-edges as a function of x 2 
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stress on the groove and tongue sides indicate that the stress 
acts in the positive or negative direction of the axis. These 
stresses are opposite to ensure the mechanical equilibrium. 
The non-symmetrical shapes and the difference between the 
maximum values on the groove and tongue sides of the 
stresses are due to the presence of the groove. 

Conclusions 

The main objective of this study was to characterize the 
interlaminar stress responsible for delamination at free- 
edges in EWF. A three-dimensional finite element model 
was used to predict the cupping deformation and to char¬ 
acterize stresses developed in the EWF. The physical and 
mechanical properties of the OSB substrate required in the 
model were determined experimentally. 

Good agreement has been found between the numerical 
and experimental EWF cupping deformation. Compared to 
the experimental measurement, the finite element solution 
overestimates the maximum deformation by 7%. Para¬ 
metric studies indicate that the maximum cupping defor¬ 
mation is influenced by the MOE in parallel and 
perpendicular directions of the OSB substrate and by the 
shrinkage coefficient in tangential direction of the sugar 
maple. The variation of these parameters with ±15% leads 
to a modification of the simulated maximum cupping 
deformation by approximately ±19%. 

The stress distribution fields generated by the model are 
in accordance with the delaminations observed on the OSB 
substrate after conditioning at 20°C and 20% RH. 
Delamination in EWF can occur principally under the 
action of the tension stress or a combination of tension and 
shear stress. Finally, simulated results show that the levels 
of interlaminar stresses are the highest near the free-edges 
of EWF. 
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